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The temporal evolution of systems undergoing symmetry breaking phase tran- 
sitions (SBTs) is of great fundamental interest not only in condensed matter 
physics, but extends from cosmology to brain function and finance |lH3l[5]. How- 
ever, the study of such transitions is often hindered by the fact that they are 
difficult to repeat, or they occur very rapidly. Here we report for the first time on 
a high-time-resolution study of the evolution of both bosonic and fermionic ex- 
citations through a second order electronic charge-ordering SBT in a condensed 
matter system. Using a new three-pulse femtosecond spectroscopy technique, 
we periodically quench our model system into the high-symmetry state, detect- 
ing hitherto unrecorded coherent aperiodic undulations of the order parameter 
(OP), critical slowing down of the collective mode, and evolution of the particle- 
hole gap appearing through the Peierls-BCS mechanism as the system evolves 
through the transition. Numerical modeling based on Ginzburg-Landau theory 
is used to reproduce the observations without free parameters. The close anal- 
ogy with other "Higgs potentials" in particle physics [4j gives new insight into 
hitherto unexplored dynamics of both single particle and collective excitations 
through a SBT. Of particular interest is the observation of spectro-temporal 
distortions caused by disturbances of the field arising from spontaneous anni- 
hilation of topological defects, similar to those discussed by the Kibble-Zurek 
cosmological model[2j. 

The behaviour of SBTs in many diverse systems is commonly studied under near- 
equihbrium (near-ergodic) conditions, where excitations on aU timescales contribute to the 
process, and the behaviour of physical quantities through the SBT is described by power laws 
and critical exponents. However, when the ordering proceeds non-ergodically, the situation 
is fundamentally different. The quasiparticles and collective boson excitations perceive the 
homogeneous crystal background as an effective vacuum. Consequently the ordering takes 
place as a well-defined sequence of events in time. For studying the temporal evolution 
of elementary excitations through a SBT in condensed matter systems undergoing second 
order SBTs, crystals with electronically-driven instabilities are perhaps the most suitable. 
Chosing a prominent example amongst such systems, rare-earth tri-tellurides have an elec- 
tronic instability caused by a Fermi surface nesting, leading to a second-order transition to a 
broken symmetry charge-density- wave (CDW) ordered state at low temperature [6]-P|. The 
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CDW state is characterized by spatial modulations ^ co8{Qx + 0) of the electronic density 
and of the lattice displacements which give rise to a complex OP ^ Ae'^^. The elementary 
bosonic collective excitations of A and are the amplitude and the phase modes (AM and 
PhM), and may be viewed as a "Higgs field opening a gap 2|A| in the fermionic spec- 
tra. Quasiparticle (QP) and AM excitations in TbTca have been recently systematically 
characterized jUni [H]. In this paper we focus on TbTe3 (Tc 336 K), but experiments 
on a number of microscopically diverse systems (DyTca, 2H-TaSe2 and K0.3M0O3) are also 
presented demonstrating the range of behavior observed within this universality class of 
SBTs. 

The main idea realised here is to repetitively quench the system into the high symmetry 
state using a short (50 fs) intense "destruction" (D) laser pulse and then monitor the time 
evolution of refiectivity oscillations using a pump-probe (P-p) sequence as it freely evolves 
through the transition (Fig. la). (The experimental details, including data processing are 
given in the supplementary information (SI)). The D pulse excites electrons and holes, 
thus suppressing the electronic susceptibility at 2/ci? whose divergence is the cause for the 
CDW formation. Any asymmetry in the band structure also leads to an imbalance of the 
e and h populations, shifting the chemical potential and causing a disturbance 5q of the 
Fermi surface kp according to — rih (x\ Sq \ /n and destroying the CDW. Following the 
quench, after initial rapid quasiparticle (QP) relaxation, we can expect the appearance of 
topologically nontrivial local configurations - domain walls, solitons. etc. which are allowed 
by the ground state degeneracy with respect to cj) [12]. Our technique allows us not only 
to control and monitor the emergence of elementary excitations and evolution of with 
high time-resolution, but permits a unique experimental detection of spatio-temporal field 
distortions arising from domain wall annihilation events. 

Focusing on TbTe3, the evolution of the system after a quench is measured by the transient 
reflectivity AR/R as a function of D — P time delay Atu (Fig. lb)). We can distinctly see 
an exponential QP transient at short times At23 < 1 ps, and an oscillatory response due to 
AM and coherent phonon oscillations [11] evolving through the SBT. A 2D plot highlighting 
the QP response is shown in Fig. 2a). We see that immediately after the quench the QP 
peak amplitude Aqp is completely suppressed, indicating the disappearance of a gap. As 
the QP peak starts to recover, initially the QP lifetime tqp is quite long (a few ps), but 
recovers quickly with time. Both Aqp and tqp recover to their equilibrium values within 
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1-2 ps. A single exponential fit to both tqp and Aqp is shown in Fig 2b), giving the QP 
gap recovery time r^Qp = ^Aqp = 650 ± 50fs. This is consistent with the previously reported 
relation tqp - 1/A(T) [HI [13]. 

Fig. Ic) shows AR/R for different At^ with the QP signal subtracted, showing an 
unusual oscillatory response through the SBT. The Fast Fourier transform (FFT) power 
spectra of these data as a function of Ati2 are plotted in Fig. 2c). The most obvious non- 
trivial observation in the uj — Ati2 plots is that the intensity of the AM fluctuates strongly 
up to ^ 7 ps. It then gradually saturates with increasing Atu (not shown). The fluctuations 
are irregular at first, showing a distinct slowing down in the critical region At 12 = 1.5 ps. 
The AM, whose frequency in the equilibrium broken symmetry state is 2.18 THz, shows a 
dramatic softening for Ati2 < 2 ps. In the course of the system recovery, the AM crosses 
the 1.75 THz phonon mode, and a Fano interference effect is clearly observed around — 1 
ps, similar to the one observed in the T-dependence.[llJ Significantly, the spectra appear 
strongly distorted around Ati2 = 3.5 — 4 ps, showing asymmetric diagonal "blobs". After 
6 ps the fluctuations die down and the AM intensity eventually reaches full amplitude in 
approximately 60 ps without any significant change in linewidth and frequency (not shown) . 

To model the evolution of the system through the SBT, we describe the field 
using a non-linear weakly dispersive Ginzburg-Landau (GL) model. Neglecting phase 
fluctuations [14j, the potential energy of the system can be described by a double- well, rather 
than a "Mexican hat" potential: 



whose time-dependence is shown schematically in Fig. 3a. Here A{t^z) = A{t^z) / A^q is 
the time-space dependent amplitude of normalized to the equilibrium value Ag^, and ^ is 
the coherence length coupling regions of different z. The function 1 — 7y(t, z) is a parameter 
describing the perturbation, akin to the temperature deviation (T — Tc) from criticality in 
usual GL theory. For spatially uniform A{t)^ r]{t) = 77(0) exp(— t/r^^^) as plotted in Fig. 
3b). Its exponential form and the parameter taqp are experimentally determined from fits 
to Fig. 2b). The time-evolution of U and fi = 1 — rj with r]{0) = 2 is shown schematically 
in Fig. 3a): Before the D pulse, and for large t or z, 77 = 0, so the system resides in a 
homogeneously ordered ground state with 1^41 = 1. Immediately after the D pulse, 1 — 77 < 
and the double well potential disappears in favor of a single energy minimum at ^4 = 0. 
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As 1 — increases and becomes positive, nonzero minima emerge at ±Amin = — vY^'^ ^ 
and start to attract the system, which is soon trapped in one of them, and the symmetry is 
broken again. From Eq. [1], the equation of motion can then be written as: 

Here Uq is the angular frequency of the bare {2kjp) phonon mode responsible for the CDW 
formation; the second term describes its damping a < /^i^am/^am- The exponentially 
decaying light intensity due to the finite penetration depth of light is accounted for by the 
excitation function rj^t^z) = r]{t) exp{—z/X) where A=20 nm is the light absorption depth 
of TbTe3 at 800 nm. Using the experimental values for rgp, uam = coq/^tt = 2.18 THz, the 
linewidth Auam = 0.2 THz, coherence length ^ = 1.2 nm [8j and penetration depth A = 20 
nm, there are no free parameters and we can compute A{t^ z). In Fig. 3b) we first plot the 
spatially homogeneous solution with <^ = 0, with and without the P pulse. In this particular 
simulation, the parameters were chosen to illustrate that a small perturbation of the P pulse 
causes the system to revert to a different minimum. With many preceding oscillations, the 
final ground state is ergodically uncorrelated with the initial one, hence the formation of 
domains is expected under inhomogeneous conditions. 

The full inhomogeneous solution A{t^ z) to Eq. 2 is plotted in Fig. 4a). We see that after 
^ 1 ps, four domains are formed parallel to the surface with A{t^ z) oscillating either around 
1 or -1 (orange or blue respectively), accompanied by the emission of A(t,z)-field waves, 
which propagate into the sample. At ^ 3 ps we observe the fusion of two domain walls, 
which is accompanied by the emission of field waves of A{t^ z) now propagating towards the 
surface and into the bulk (arrows). They appear to reach the surface around Ati2 — 4-5 
ps which - as we shall see - cause detectable distortions of the spectra at around 5-6 ps. 
The appearance of such 74-waves following annihilation is qualitatively robust with respect 
to the parameter values for TbTcs (see SI). (More 74-field wave dynamics is shown in the 
accompanying movies.) 

The calculated refiectivity response detected by the probe p is given by the dif- 
ference between the response with and without the pump P pulse: AR{t^ Ati2) oc 
f^[A'j^{t^ z) — A'j^p{t^ z^ Ati2)]e~^/^dz^ where A'j^p{t^ z^ Atu) is calculated replacing rj 
Vp{t) = '^(0) exp(— t/r^p) + HB(t — ti2)exp[—{t — ti2)]/rsp). The exponential term accounts 
for the probe penetration depth. The typical value of H = 0.1, where H is the P pulse 
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intensity relative to the D pulse and B(t — ^12) is a unit step function (see SI for details). 
The calculated response for the homogeneous solution Ai?(t, Ati2) is shown by the green 
curve in Fig. 3b). 

In Fig. 4b) we show the FFT power spectra from AR{t^ At^) taking full account of 
spatial inhomogeneity for D, P and p. The main features of our data in Fig. 2c) are 
unmistakably present: oscillations of A{t) are clearly visible at short times, as well as the 
hallmark of the transition itself, namely the critical slowing of the AM oscillations close to 
the critical point tc — 1-5 ps, pinpointing the exact critical time of the transition tc- At 
this bifurcation point topological defects are formed. The calculation also reproduces the 
softening of the AM for Ati2 < 2 ps. After 2ps, the ripples in A{t^ z) discussed above cause 
a temporal deformation of the spectral profiles, giving diagonal blobs at 5 ^ 6 ps shown in 
Fig. 4b). These are remarkably similar to the diagonal spectral distortions observed in the 
experimental data in Fig. 2c). A summary of exhaustive modeling within a wide parameter 
space is presented in the SI showing that inhomogeneity without A-field waves cannot cause 
diagonal deformations, only vertical ones. The diagonal distortions muj — tu plots are thus 
unambiguously attributed to the A-field waves created upon the annihilation of defects. 

Identical experiments on three additional microscopically diverse systems (^i7-TaSe2, 
K0.3M0O3 and DyTe3) displaying a 2nd order SBT presented in the SI show that the se- 
quence of events after the quench: (i) ultrafast QP gap recovery (ii) -field amplitude 
fluctuations (Hi) critical slowing down through tg and (iv) domain creation coherent 
defect annihilation is commonly observed in systems which unambiguously belong to the 
same universality class (the tellurides and the selenide). The microscopic properties of the 
underlying vacuum such as A and ^ change the details. K0.3M0O3 - which already shows 
lack of AM softening in the T-induced transition[13j- does not show step (iv), which we 
attribute to departure from universality. It is interesting to note that the mechanism de- 
scribed here for topological defect creation is conceptually and historically related not only 
to vortex formation in superconductivity, but also to the Kibble-Zurek mechanism for the 
formation of cosmic strings. The A{t^ z) -waves such which we observe after annihilation 
events have a direct analogue in the Higgs spontaneous symmetry breaking mechanism. All 
these models share a common underlying potential, albeit with different symmetries of OP 
and microscopic properties of the underlying vacuumfH EJ HI [T5j. A notable distinction of 
our system is that (j) relaxation is slow compared to the relaxation of the potential itself. 
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allowing the collective mode and topological defect dynamics to be clearly observed. In 
superconductors, the QP relaxation is beautifully observed p!6]. but the collective mode is 
overdamped and unobservablepT] . Thus real-time observations and coherent control of the 
macroscopic order through the SBT appears to be a rather unique feature of femtosecond 
laser experiments on CDW systems. 
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Figure 1: a) A schematic diagram of the timing of the laser pulses: a destruction (D) pulse 
(represented by the blue ball) quenches the system, while a pump-probe (P-p) sequence probes the 
reflectivity at a later time Ati2. P and p pulses are represented by red and green balls respectively, 
b) Raw transient reflectivity data AR/R for different delays At 12 (displaced vertically), showing a 
QP peak at short times, and OP and coherent phonon oscillations at longer times c) AR/R with 
the QP response subtracted. 
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Time after quench Atn [ps] 

Figure 2: \/\R/ R\qp after the quench as a function of At 12, (Note the ripples arising from coherent 
oscihations of the OP.) b) The QP hfetime tqp and the amphtude of the QP response Agp as 
a function of Ati2. A single exponential fit to both data sets (shown by the lines) gives TrQp — 
^Aqp — 650=b50fs. c) The FFT power spectra of the data in Fig. 1 c) as a function of Ati2 recorded 
at 100 fs intervals. Note the non-periodic fluctuations of intensity at around the transition (1.5 ps) 
and the strongly asymmetric space-time lineshapes for Ati2 = 2 ^ 4 ps (white arrow). The orange 
arrow indicates the critical time of the SBT. The red line is a superimposed plot of a fit to the QP 
decay from b) above. 
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t < t ~ 50 fs t ~ 650 fs t I ps t > 6 ps 




Figure 3: a) The evolution of the potential U (Eq.l) as a function of time. The system is in the 
high symmetry state (grey) at very short times after the quench. The red dot signifies the state of 
the system, while the blue/orange potential signifies a topologically mixed broken-symmetry state, 
b) The time-dependence of the control parameters = 1 — r]{t) (solid red) and /ip(t) (dashed 
red) calculated with the experimental value of ta^p = 0.65 ±0.05 ps forAti2=0.4 ps. The predicted 
oscillations of A(t) with and without the P pulse are shown by the dashed and solid oscillatory blue 
curves respectively. The predicted optical response AR(t) is shown by the green curve. 
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Figure 4: a) The calculated A{z^ t) as a function of depth z and Ati2. Note the ripples in caused by 
the annihilation event at ^3.5 ps (arrows), b) The corresponding computed transient reflectivity 
AR(z,t) as a function of Ati2. Note the predicted diagonal distortion due to the A{z,t)-wdive 
reaching the surface, indicated by the white arrow. The orange arrow points to the critical slowing 
down at the critical time of the transition tc, i.e. the bifurcation point. The thin red line shows the 
function r]{t). 
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